The route of frustrated cuprates from antiferromagnetic to ferromagnetic spin-1/2 
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From thermodynamics, LSDA+17 studies and exact diagonalizations of a five-band Hubbard model 
on Cu02 stripes we found that L^ZrCuCU (Li2 0uZrO4 in traditional notation) is close to a ferromag- 
netic critical point. Analyzing its susceptibility x{T) and specific heat c p (T, H) within a Heisenberg 
model, we show that the ratio of the 2nd to the 1st neighbor exchange integrals q=— J2/Ji~0.3 is 
close to the critical value a c =l/4. Comparing with related chain cuprates we explain the rather 
strong field dependence of c p , the monotonous down shift of the peak of x{T), and its increase for 
a— >a c +0. 



The one-dimensional (ID) spin-1/2 antiferromagnetic 
(AFM) Heisenberg model (HM) is one of the most stud- 
ied many-body models in theoretical physics. Much of 
its physics is now well understood based on the rigorous 
Bethe-Ansatz method for infinite chains [l[ and on finite 
cluster calculations. Thermodynamic benchmarks of this 
model relevant here are: (i) single maxima of the spin 
susceptibility x(T) at fc B T m w 0.64J and of the specific 
heat c v (T) at fc B T m w 0.48 J, (ii) c„ cx T/J at T -> 0, and 
(hi) X *(0) = J x (0)/Ng 2 [il = 1/tt 2 and d X (T)/dT -++oc 
at T — > 0. Hereafter J = J\ denotes the nearest neigh- 
bor (NN) exchange. For ferromagnetic (FM) J\ < 0, 
X(T) oc 1/T 2 and c v oc y/T/ | Ji | at T -> 0; c v shows 
a broad maximum at k B T^ = 0.35 | J\ | and a field in- 
duced 2nd maximum at low T and H < 0.008 | J\ \ /gfJ- B 
■ The general Hamiltonian 7i with next-nearest neigh- 
bors (NNN) J2 or further in-chain exchange Ji included 
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has also attracted attention due to the frustration caused 
by AFM J2, irrespective of the sign of J±. If the Ji 
are AFM, the frustration may cause a spin gap, e.g. 
for J2/J1 > 0.241 and Ji = 0, i > 3 (adopted mostly 
below). It strongly supports a dimerized ground state 
in spin-Peierls chains such as in GeCu0 3 1- Re- 
cently, FM-AFM analogs realized in most edge-shared 
chain cuprates have caused attention with respect to 
strong quantum effects to unusual thermodynam- 
ics of the disordered phase d, 0, 01, and to helicoidal 
ground states found in some chain cuprates at low T 
E B E El, O, El 0, El, E, E3 ■ However, issues like 



the behavior at very low T and in magnetic fields near the 
critical point a c = —J2/J1 = 1/4 arc still unclear and dif- 



ficult to study numerically [fj] even by the transfer matrix 
renormalization group (TMRG) method. For a > a c the 
ground state of a classical chain is formed by a helix with 
a pitch angle <fi given by cos0 = — J1/AJ2 = 1/(4q). This 
helix interpolates between a FM-chain at < a < a c and 
two decoupled AFM-chains at a = 00. Noteworthy, a c 
is unaffected by quantum effects 17|. Since this should 



hold for the case of long-range inchain couplings, too, we 
expect a down(up)shift of a c for AFM (FM) J,, (i > 3): 
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Recently, low T- X (T) data for Rb(Cs) 2 Mo 3 Cu 2 Oi2 [15, 
have been refitted by the isotropic J1-J2 HM near 
a c . However, both compounds seem to be affected by 
Dzyaloshinskii-Moriya interactions Djj (Si x S 3 ) and 
exhibit a very complex, partially unresolved crystal struc- 
ture, complicating a theoretical study even more. 

Hence, studies of less complex systems described by 
Eq. (1) but with | a — a c |<C 1 are of general interest. 
Analyzing x(T), cp(T,H), and the electronic structure 
of Li2ZrCuC>4 we will show that it is a suitable candidate 
to probe the vicinity of a c from the helical side. Together 
with data for related systems with a > 1 it provides a so 
far missing link near a c to study e.g. the a-dependencc 
of relations (i)-(iii), moving from AFM to FM-chains. 

The orthorhombic crystal structure of Li2ZrCu04 
[IH (space group Cccm) with the lattice constants 
a=9.385 A, b=5.895 A, c=5.863 A is shown in Fig. Q] 
Here chains (formed by flat edge-shared CUO4 tetrahedra 
like the edge-sharing of CUO4 plaquettes in other chain 
cuprates) run along the c-axis. Also the Cu-0 bond 
length of 2.002 A and the Cu-O-Cu bond angle 7 = 94° 
resemble those with FM Ji. 
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FIG. 1: (Color online) Crystal and electronic structure near 
the Fermi level -E F =0 of Li2ZrCu04. Left: Crystal struc- 
ture, Cu 2+ - large orange o, Zr 4+ - light magenta o inside 
the magenta corner-shared ZrOe octahedra, red o- O 2- , and 
light blue o- Li + (Li (split) positions near Zr are omitted 
for clarity); the nonplanar edge-shared CUO2 chains (olive- 
green). Right: LDA-FPLO band structure (o) and TB fit (red 
line). F,X,Y,Z,S are symmetry points in wave vector nota- 
tion: (0,0,0);(27r/a,0,0);(0,27r/6,0); (0,0,27r/c);(27r/a,27r/b,0), 
respectively. 



The 7-polymorph of Li2ZrCu04 (Li ordered) was pre- 
pared by a solid state reaction of L12CO3, ZrC-2 and CuO 
U9fl. The reagents were mixed in an agate mortar and 
fired for a few hours in a Pt boat at 700°C to decarbonate 
them. Final firing of the pellet was performed at 1050°C 
for 24 h in a flow of O2 followed by furnace cooling in 
O2. Phase purity was confirmed by x ray diffraction. 

The magnetization of Li2ZrCu04 measured in a range 
2< T < 350 K for 0.1 T by a Quantum Design SQUID 
magnetometer is shown in Fig. 2. From the observed 
T m w 7.6 K one might at first glance expect an AFM 
spin liquid regime with J\ or J2 ~ 12 K, if 7 is just 
by chance close to that bond angle where J\ changes 
its sign and either J\ ;2> J2 > due to the non- 
ideal chain geometry or vice versa J2 3>| J\ |. But 
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FIG. 2: (Color online). Magnetic susceptibility of Li2ZrCu04 
together with fits by the J1-J2 model for periodic chains. 
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FIG. 3: (Color online) Specific heat c p of Li2ZrCu04 vs. T at 
various external magnetic fields H. Inset: the same for c v /T. 



the measured x*(T m ) is twice as large as the AFM- 
HM value of 0.1469 (x AFMHM (T m )=0.0183 emu/mole for 
g=2). l/x(T) oc T + Q C w reveals a FM Curie- Weiss 
temperature Qcw= -24 K using a narrow temperature 
range near 350 K. Both facts exclude any AFM-HM like 
scenario. But they point to FM exchange involved in 
accord with fits by the J1-J2 model (Fig. 2). 

Specific heat down to 0.35 K was measured by Quan- 
tum Design Physical Properties Measurements System 
(see Fig. 3). It shows a relative sharp peak near 6.4 
K at H =0. Using c p ~ c v = c [2fJ. the observed 
ratio T n j/T^ 1 = 1.17 differs from 1.33 predicted by the 
AFM-HM. Note that T m nearly coincides with the T 
for which d%(T)/dT becomes maximal. Hence, it is 
unclear whether can this peak be attributed either to 
a c p anomaly indicating often a magnetic phase transi- 
tion [2l| , or to a specific feature of the disordered phase 
generic for the ID frustrated J1-J2 HM at ctc< en < 0.4 
. Here c v exhibits a iwo-peak structure @, a sharp 

peak at low T under consideration and a broad one at 
high T hidden in the phonon region ( k B T ~ 0.65 | J\ \~ 
260 K in the present case) . Anyhow, with increasing field 
T m is downshifted and c p {T m ) is suppressed but c p {T) in- 
creases rapidly for T > 12K, well above a possible phase 
transition near 6K. 

A similar strong iJ-dcpcndence is found in full diag- 
onalization studies of large rings, where the low-T peak 
is first downshifted with increasing H and upshifted at 
higher H (Fig. 4). The strong iJ-depcndcncies of both 
Ac(H)=c(H)-c(0) and -AT m = T m (0)-T m (#) oc H 2 al- 
ready at weak fields H < 9T result from the vicinity to 
ct c [22j. Adding a usual lattice contribution c\ a t oc T 3 
(G>£)=220 K) to the calculated spin specific heat within 
the isotropic J1-J2 HM the experimental data are best 
described by cv = 0.3 (Fig. 4). From the low ff-crossing 
point near 12 K we estimate J\ « 405 K to ~ 363 K 
using the x(T) data for a = 0.29. The low-T peaks ex- 
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FIG. 4: (Color online) Specific heat c„ vs. temperature in 
units of | Ji | for various magnetic fields H within the FM- 
AFM-J1-J2 HM for a = -J 2 /Ji=0.3 and a ring with N=20 
sites (upper row). Finite size effects for the maximum of x(T) 
(left) and the low-T maximum of c„(T) (right) (lower row). 
The N = 00 values are taken from TMRG results of Ref. H 



trapolated to N = 00 would be expected near k^T^ « 
0.013(0.0115) I Ji I, respectively, i.e. near 5.3(4.2) K be- 
low the observed one at 6.4 K, similarly as the expected 
T m « 4.5± 1.7 K is below the observed one near 7.6 K 
(Figs. 2-4). These observations are in accord with the 
scenario of a phase transition at 6 K as discussed above. 
The slightly different results from fitting x(T) or c p (H, T) 
might be due to anisotropics and interchain coupling. 

To estimate the interchain coupling, wc consider the 
measured Curie-Weiss temperature Qcw — rQcw= - 
24 K, where 1 > r(T) « 0.25 is estimated from the 
calculated dx~ 1 (T)/dT taken at the highest available 
T. Here it is still outside the asymptotical CW-range 
fc B T Ji , where r — ► 1. The high-T expansion of 
X (T) yields S C w = 0.25 i. e. 



28c.iv/r = Ji (1 - a) + J 3 + J ± + 2 J dl + 2 J d2 



(3) 



where the neighbor number Zi = 2 for couplings along the 
c and 6- axes and z,;=4 for diagonal interchain exchange 
(dl,d2) within the 6, c-plane. From the tight binding 
(TB) fit of the band dispersion we find similar direct 
and diagonal interchain transfer integrals t± ~ tdi ~ t^i ■ 
Setting Ji = Jdi = Jd2, we found | J\ , J2 S> J± ~ 9 



K in accord with the LDA results (J±~ 7 K). Thus, the 
adopted ID magnetic approach is a reasonable starting 
point despite the more 2D electronic structure seen e.g. 
along the symmetry lines T-Y and S-X in Fig. 1. 

To get insight into the J-sct obtained above, we 
performed calculations of the electronic and magnetic 
structure within the local (spin) density approximation 
(L(S)DA). In addition, LSDA+t/ calculations and exact 
diagonalizations for an appropriate extended multiband 
Hubbard model were carried out to take the strong corre- 
lation for the Cu 3d holes into account. The LDA calcula- 
tions (Perdew-Wang92 parametrization) were performed 
using the full-potential local-orbital minimum-basis 
scheme (FPLO, vers. 5.00-19) [2^. We employed a basis 
set of Cu(3s3p):(4s4p3d), 0(2s2p3d), Zr(4s4p):(5s5p4d), 
and Un S ):(2s2p3d). For the LSDA+C/ in the AFM ver- 
sion [24j we used E/3<j=6.5±1.5 cV and the intraatomic 
exchange J ex =l eV. Following the approach of analyzing 
total energy differences for various magnetic superstruc- 
tures we obtain J x = -151^35 K and J 2 =35t12 K. 
Using a typical one-band Hubbard U e jf ~ 3.5 eV as well 
as t<z and t$ from the TB fit of the band at E F (Fig. 1) 
results in J2=46 K and J 3 =6 K employing J; = Atf/U. 
Thus, we arrive close to cv c =0.195 in the present case of 
J 3 ? (s. Eq. (2)). 

Finally, a collection of known T^/ J 2 and x(^nr) val- 
ues from other chain cuprates we derived from their x{T) 
data @, 0, EH H Hjl, is shown in Fi g- 5 US- In 
particular, it is clear why the large-a chains in SrCu02 
and LiVCu04 are often regarded as AFM-HM archetypes 
[jjl ]. Only after the discovery of spirals, detailed inelas- 
tic neutron scattering studies, and our three component 



a. = 00 (J, = 0), 1D-AFM Bethe-Ansatz 
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FIG. 5: (Color online) Empirical T m in units of the fitted 
J2 value of the FM-AFM Ji-J2-model for several frustrated 
chain cuprates (black squares): 0- a c =l/4, 1- Li2ZrCu04, 2- 
Pb2[CuS0 4 (OH) 2 ], 3- RD2CU2M03O12, 4- Cs 2 Cu 2 Mo 3 Oi2, 5- 
LiCu 2 2 , 6- NaCu 2 02, 7- LiVCu0 4 , and 8- SrCu0 2 . The 
measured pitch is given in brackets. The filled green o denote 
the full diagonalization results of the J1-J2 model on rings 
with JV=20 sites. Inset: the maximum value of x(T). 
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theoretical analysis (HM, Cu-0 Hubbard model, LDA) 
initial assignments for L1VCUO4 and LiCu202 were cor- 
rected d, [3, E3j E3 ■ Similarly, among systems assigned 
so far as "perfect" realizations of the AFM/FM HM (e.g. 
[UHil) could be further J1-J2 candidates worth to be re- 
visited. Similar plots which accent the FM critical point 
can be made also for the low-T peak of c v or %(0) (which 
monotonously increases and diverges finally as a — * a c ) 
[30j | . We expect that a vanishing Tf^ and a diverging 
%(T) for T — > in approaching a c are generic for a FM 
critical point. It should hold also for models beyond the 
J1-J2 HM. Further couplings do affect the helical phase 
in changing e.g. the pitch and a c (see Eq. (2)). 

Comparing a ~ a c for I^ZrCuC^ with a > a c we 
found for other chain cuprates the question arises, what 
is the microscopic reason for? There are at least two 
options: (i) an enhanced | J\ | at a standard J2 value 
and vice versa (ii) a slightly enhanced | J\ | at a re- 
duced J 2- Case (i) can be ascribed to enhanced FM 
contributions to J\ which arise from the direct exchange 
K p d or from the Hund's rule coupling at the sharing 
O ions within a five-band Cu 3c? O 2p extended Hub- 
bard model. Unfortunately, there is no generally ac- 
cepted K p( i value; but it is the most sensitive qua ntity 
for the determination of Ji in edge-shared cuprates [3l( . 
Nevertheless, usually K p d is treated as a fit parameter: 
The well-studied Li2Cu02 can be described with K p d=50 
meV [Hj], whereas microscopic calculations for La2Cu04 
yield 180 meV [33| and a structural analysis of GeCu03 
was performed adopting K p d=ll0 meV [31j. Within the 



Cu 2>d yz O 2p y ,p z extended Hubbard model for planar 
Cu n 02n+i open chains (n < 5), we adopted K p d=70 
meV. From a direct mapping onto the J1-J2 HM using 
two O site energies A Py d=2.7eV and A Pi d= 3.2eV as well 
as Li 2 Cu02-like parameters, we found Ji=-317 K, J 2 ~ 
90 K, and a =0.284, close to our empirical values. In 
case (ii) supported by the LSDA+L 7 results, we arrive 
also close to a c - Here J2 amounts 46 K, only. From 
a comparison with other cuprates in Fig. 5 more insight 
will be gained into the nature of the exchange in cuprates 
and into the novel physics generated by a quantum FM 
critical point. 

To conclude, we have shown that a growing number of 
edge-shared chain cuprates form a special family which 
thermodynamics can described within the J1-J2 model 
with FM NN and AFM NNN exchange. Moving from the 
AFM-HM towards a c , almost achieved for I^ZrCuCU, 
observed monotonic changes can be explained. Only 
chains near the FM critical point show peculiar phys- 
ical properties such as the strong iJ-dependence of c p 
in a large T range reported here. Further studies of 
Li2ZrCu04 at very low T, under pressure, and in high 
fields are highly desirable. If the observed c p -peak is re- 
lated to magnetic ordering, neutron diffraction below 6 
K should reveal a spiral with a pitch below the minimum 
value of 62° observed so far among edge-shared chain 



cuprates for LiCu2 02 [8j. Inelastic neutron scattering 
studies might be helpful to refine the exchange integrals, 
especially with respect to the interchain coupling. 
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